function Fig2_d()
lb = 25.66 / 6^0.5;
EC = 56.1 * 6^0.5;

% This proccessed data is processed by 
%  is located at .../processed data
result = readmatrix('processed data/result_one_half.csv');
r = result(:,2) * lb;
phi_r_total = result(:,3) * EC;
phi_r_imp = result(:,4) * EC;
phi_r_2DEG = result(:,5) * EC;
n_r_2DEG = result(:,6) / lb^2;


Fcolor='k';
Bcolor='w';
Fsize=18;  %font size

f = figure();
hold on
plot(r,phi_r_total)
plot(r,phi_r_imp)
plot(r,phi_r_2DEG*-1)

f.Color = Bcolor;
ax = gca;
box on
ax.Color=Bcolor;
ax.XColor=Fcolor;
ax.YColor=Fcolor;
xl=xlabel('r (nm)');
yl=ylabel('(meV)');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;
axis([0,125,-18,2])

end
